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We study a generic counting process for systems governed by ageing waiting times. In contrast 
to renewal continuous time random walks with independent waiting times, we consider the case 
when each state is characterized by its own internal waiting time process, having been initiated at 
time t = 0. Therefore each transition from state n to n + 1 is triggered by ipi, the forward waiting 
time density, instead of the regular waiting time density ip. For states characterized by heavy-tailed 
forms of ij) we obtain an asymptotically logarithmic evolution of the counting dynamics n(t), whose 
fluctuations vanish relatively to the mean. The counting process introduced here describes the 
dynamics of crack propagation, the counting of on-off transitions in ensembles of blinking quantum 
dots, or, more generally, the dynamics in complex systems such as glasses or biological cells. 

PACS numbers: 82.20.-w, 87.10. Mn, 02.50.-r, 05.40.-a 



In the classical Scher-Montroll- Weiss continuous time 
random walk (CTRW) model a new transition in the sys- 
tem, for instance, a nearest-neighbor jump on a lattice, 
occurs after a random waiting time [lj, |2|. Individual 
waiting times are independent and drawn (renewed) from 
the common probability density tp( T )- Power-law forms 
V>(t) ~ T -i-a W1 t n o < a < 1 were originally intro- 
duced for charge transport in amorphous semiconductors 
2]. In absence of a bias tp(r) gives rise to subdiffusion 
(x 2 (t)) ~ K a t a , where K a is the anomalous diffusion co- 
efficient d,|3|- Such long-tailed forms of ip(r) were shown 
to govern the motion of tracer particles in the cytoplasm 
Q and in membranes @ of living cells, in reconstituted 
actin networks Q , as well as the blinking dynamics of sin- 
gle quantum dots [1| and the dynamics involved in laser 
cooling 11]. Physically, the form ip(r) may arise from 
comb models [§] or random energy landscapes [ToJ • The 
divergence of the characteristic waiting time J Q tiP(t)oIt 
leads to phenomena such as ageing [l2j and weak ergod- 
icity breaking [13}, with profound consequences for our 



conceptions, e.g., in molecular cellular processes [14 1. 

Here we relax the renewal property of CTRW pro- 
cesses and consider the counting process n(t) in a system 
in which waiting times are no longer a property of the 
walker in an annealed environment, but at the outset a 
property of each state n. The random quantity n could 
be the number of jumps in between independent CTRW 
processes in a completely biased random walk, see Fig. Q] 
At each n, events that can trigger a transition occur ac- 
cording to independent waiting time processes governed 
by ip( T ) — with arbitrary exponent a, and we as- 

sume that all these processes were initiated globally at 



time t = 0. As triggering events occur independently at 
different states, we typically arrive at n in between two 
triggering events. Arriving at n at some time t' , the time 
span t until the transition to n + 1 is triggered is thus 
given by the forward waiting time density iJ)i(t\1/), which 
explicitly depends on the ageing time t' [l2|. In the the- 
ory of renewal CTRW, f/'i( T l*') governs the occurrence of 
the first jump in an aged system [12, EE [Hf, while all 
successive waiting times are drawn independently from 
"0(t). We here consider a process in which each transi- 
tion is controlled by xpi(r\t') and quantify the probability 
h n (t), and its associated moments, to find the system in 
state n at time t. Interestingly, we find a logarithmic 
growth of the moments with time and that the distribu- 
tion h n (t) becomes increasingly sharp for < a < 1 and 
a > 2, while its width remains comparable to its mean 
in the intermediate regime 1 < a < 2. 

As an example, consider the propagation of a linear 
crack front in a solid. At each point n measured in the 
propagation direction of the crack, a vacancy is moving 
diffusively perpendicular to the n-axis. If the crack has 
propagated up to site n — 1 and the vacancy at n reaches 
the n-axis, the crack grows to n, etc. The local triggering 
process is thus given by a waiting time density of the form 
ip(r) with exponent a = 1/2 corresponding to the return 
of a one-dimensional random walk At increasing n, 
and thus increasing time, the local triggering process on 
average ventures further away from the n-axis, such that 
waiting times typically increase with n: the process ages. 
In a complex system the triggering vacancy could itself 
diffuse with a power-law waiting time density with expo- 
nent < 7 < 1, such that a — j/2. Alternatively, in an 
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FIG. 1: Ageing counting process: propagation among suc- 
cessive states n as governed by the transition time density 
p n (t). At each n triggering events occur according to the 
waiting time density i[>(t), globally initiated at t = 0. After 
the walker arrives to n at time t' , it has to wait for the next 
triggering event to be allowed to continue to state n + 1 after 
the forward waiting time r. 



ensemble of quantum dots or tracer particles in an actin 
network, we could be interested to count the number of 
events when the next quantum dot turns from an on to an 
off state or when the next tracer particle breaks out from 
its cage. As final examples we could consider the dynam- 
ics of train hopping in England, where the delay times 
indeed follow power- laws with exponent a ~ 1.8 fl7j j . 
or message relaying along human communication paths. 
For instance, the Einstein and Darwin correspondences 
had delay times characterized by a s» 0.5 [18[. A simi- 
lar scenario could apply to cascades of active molecular 
transport in biological cells. The non-renewal counting 
process considered here is a general tool to quantify tran- 
sitions in complex systems, in which such events depend 
on local triggering dynamics governed by power-law dis- 
tributions of index a > 0. 

To formulate the model mathematically we introduce 
the probability density p n (t) to arrive at state n at time 
t. This quantity fulfills the recursion relation 



Pn(t) 
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Jo 



P n-l{t')Mt-t'\t')dt', 



(1) 



where ip\ (r\t') is the density of the forward waiting time r 
that the walker has to spend in a new state after having 
arrived there at time t'. If the local triggering process 
follows ip(r) with < a < 1 (see below for the case a > 
1) then the waiting time r from time t' until transition 
to the next state is distributed according to the density 



sin(7ra) 



T a (t' + t) 



(2) 



for the first jump in ageing CTRWs 0, [H, E0] ■ If 
the walker is initiated in state n = at time t = to, the 
recursion (TTJ is closed by po(t) = S(t — t ). With respect 
to a short waiting time cut-off tq in the power-law form 
of ip(r), Eq. @ is the asymptotic result for t' tq. 



Introducing the function A(x) = i/ji(x — l|l)0(x — 1) 
with the Heaviside step function 9(x), Eq. (fTJ) becomes 



Pn(t) 
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(3) 



This expression can be interpreted as the transforma- 
tion of the product of two independent random vari- 
ables: denoting by t n the arrival time at state n, Eq. (|3]) 
states that t n = Xnt n -i, where Xn is independent of 
t n -i and has distribution A(x). In retrospect this is 
not surprising since t n _i is the only time scale deter- 
mining t n [see Eq. (|2|)]. This means that we may write 

in = Xntn-i = Xn'-'Xlh, Or \n(t n /t ) = J2i=l^ n Xi- 

Thus, in logarithmic time the transitions have indepen- 
dent waiting times. At large n we thus apply the central 
limit theorem. For 2, = In Xi we have the distribution 



(S(z -h)) = e z X(e z ) = 



sin(7ra) 9(z) 
7T (e z - l) a ' 



(4) 



with the moment generating function G(p) = (e pZi ) = 
T(a — p)/[r(a)r(l — p)]. This defines the cumulants 
K m = {d/dp) m lnG(p)| p =o- In particular, we find the 
average p, = k\ = —T'(a)/T(a) — 7 and the variance 
cr 2 = k 2 = -tt 2 /6 + d 2 lnT{a)/da 2 , where 7 = ~r'(l) is 
Euler's constant. For large n wc then infer from the cen- 
tral limit theorem that ln(i„/io) is normally distributed 
with average np and variance na 2 . Thus, for large n the 
approximate distribution assumes the Gaussian form 
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(5) 



To proceed more systematically we apply a Mcllin 
transform f{p) = J °° ff- 1 f(t)dt to Eq. ©. By the 
Mellin convolution theorem, Eq. ([3]) leads to p n (p) — 
X(p)p n -i{p) in terms of the Mellin variable p, such that 
A(p) = G(p — 1). With po(p) = ^o _1 we obtain p n (p) — 
Hp) n t P ■ While no simple expression for the inverse 
p n {t) exists, it can be shown that the asymptotic behav- 
ior at t — ¥ 00 is given by 
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(6) 



Thus, the distribution of arrival times at state n mirrors 
the power-law tail of the original distribution ip(r) plus a 
logarithmic correction with power n. Unsurprisingly the 
associated characteristic time scale diverges. 

We are now ready to turn to the central quantity in 
this study, the probability distribution h n (t) to find the 
system in state n at time t, equivalent to having arrived 
at n at t' < t, and not having moved since, 



h n (t) 



Pn (t') / MT\t')dT'dt'. 



(7) 
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FIG. 2: Average counting process (n(t)) versus time t. Sym- 
bols represent simulations for various a, the dashed lines show 
the asymptotic behavior, Eq. (|10[) for q = 1. Simulations used 
to — 10 3 and To = 1, with averages over 10 7 ensembles. Inset: 
convergence to the theoretical results with to for a — 0.75. 



The average state number (n(t)) at time t and cor- 
responding higher moments thus become (n q (t)) — 



' J h n (t). The Mcllin transform of Eq. Q is 



hfiip) = Pnip + - l]/p, and therefore h n (p) = 

t^ ) G(p) n [G(p) — l]/p. To find the asymptotic behavior of 
the moments we note that G(p) is an increasing function 
of p (as long as p < a) which grows larger than unity 
when p is larger than zero. Thus the moments (n q (p)) 
diverge as p — > 0~, and therefore have the fundamental 
strip — oo < p < 0. The long time asymptotic behavior 
of the moments is therefore dominated by the singularity 
at p ~ 0. Expanding G(p) at p ~ to second order, 
G(p) ~ 1 + fip+ \ {<J 2 + fj, 2 )p 2 . If we similarly expand the 
sum over n in (n q (p)} in deviations of G(p) from unity, 
making sure to keep the lowest and next to lowest order 
contributions, wc find 



G(p) q q\ 



(q - l)G(p)«- 



[l-G(p)]^ 1 2[1 -G(p)]« 



(8) 



Collecting terms we finally obtain the expression 

r(« + i)*g 



(n'(p)> 



/z?(-p) 



+ | [m-^/m]}, p^o- 



(9) 



Inverting the Mellin transform, we find that at large t 



(n q (t)) 
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1 



In(t/*o) 



(10) 

The first two moments to leading order grow like (n(t)) ~ 
ln(t/to)//i and (n 2 (t)} ~ ln 2 (i/A )/V 2 ! respectively. 

Fig.[2]compares the mean number (n(t)} (fT0| with sim- 
ulations, for the concrete form - 0( T ) = ar oV( T + T o) 1+Q - 
Except for a — > 1 the simulations agree excellently with 



FIG. 3: Standard deviation versus mean as a function of time. 
Simulations (symbols) are compared to the asymptotic results 
for large t [Eq. (|11[) . dashed lines]. Parameters used: to = 10 3 , 
To = 1, averaged over 7 x 10 5 ensembles. 



the asymptotic result (|10[) . From the inset of Fig. [2] we see 
that the mismatch is due to the fact that to is not large 
enough, i.e., compared to to the distribution ^>i(r|to) has 
not reached its asymptotic form, Eq. 

In Eq. (|10p the g-dependence of the dominating term 
for the moments appears naturally, corresponding to a 5- 
function behavior for the limiting distribution. From the 
subdominant terms of (n q (t)} we see that the deviation 
divided by the mean decays logarithmically as 



y/(n*(t))-(n(t))< 
<»(*)> 



(11) 



as nicely corroborated from simulations in Fig. [3l This 
implies that the standard deviation of the counting pro- 
cess versus the mean number (n(t)) becomes increasingly 
sharp, in contrast to the behavior of the position coor- 
dinate in biased subdiffusive CTRW processes, in which 
the ratio of deviation and mean tends to a constant Q ■ 
To obtain the full form of the distribution h n (t) we 
note that its Mellin transform h n (p) can now be written 
as h n {p) = 1? Q e nXnG{ P)[G(p) - I}/ p. In the large n limit, 
h n is different from zero only for p R3 0. To obtain an 
approximate result for h n (t) we expand h n {p) for large n 
and small p, keeping the product np v constant, where the 
scaling exponent v is chosen as small as possible while 
still obtaining a non-trivial result when discarding the 
small terms. For v = 1 a S- function is obtained for h n (t), 
for v — 2 we find to zeroth order in small quantities 



hn \p) = A^o exp[n([ip + a 2 p 2 /2)}, thus 



V27 



: exp 
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2<r 2 n 



(12) 



This expression systematically improves by inclusion of 
higher order terms in p and np 3 . To first order, h!n (p) — 



4 



0.09 



0.06 



0.03 



Simulation 



0.25 







15 



Q ab6oc 



30 



FIG. 4: Probability of being in state n at a given time t. 
Lines: analytical results from Eqs. (|12|l and (|13p . circles: 
results from simulations, averaged over 7 x 10 ensembles. 



Parameters used: t — 7.9 x 10 
a = 0.25. 

h^(p)[l + (^+^)p/(2^ 

r 2 



to 



r = 1, and 



«3np /6], thus 




K 3 n 



-(y 3 -3y) 



6 ( CT 2 n )3/2 ' 

(13) 

In Fig. |4] we compare 



where y = \\n(t/ta) — fJ.n]/y 

hn\t) with simulations, demonstrating very good agree- 
ment for the dominating part of the distribution. 

What is the behavior when a > 11 In this case 
ipi(r\t') has a finite limit ipi (r) = limf-^oo V , i( T l*') 
given by ipi( T ) = ( r ) 1 4>( T ')dT' [2l| . where (r) = 
Jq 00 Tip(r)dT. Thus at long times one obtains "01 ( T ) ~ 
(q - l)A/[(r)r Q ]. For 1 < a < 2 this means that the 
counting process corresponds to the propagation as a one- 
sided Levy flight on the time axis. Inverting this to find 
the distribution of the number of transitions to reach a 
given time will correspond to solving a first passage prob- 
lem for one-sided Levy flights. From this argument we 
find that the moments are given by 19|, 22 1 



(n«(i)> 



r(g + l) / (r)r- 



T(q(a-1) + 1) \T(2-a)A 



(14) 



Thus for 1 < a < 2 the fluctuations grow as fast 
as the mean n. For a > 2 we replace a — > 2 and 
r(2 - a)4 -> (r 2 )/2, such that (n«(<)> - (2(r)t/(r 2 ))' z . 
In this case the deviations vanish relatively to the mean. 
Thus, the deviations become increasingly negligible com- 
pared to the mean value (n(t)) when the exponent a is 
smaller than unity or larger than 2. In contrast their 
ratio saturates in the intermediate range 1 < a < 2. 

The results we presented in Eqs. (O to (fTB"]) apply more 
generally, beyond the ageing waiting time process con- 
sidered so far. In fact, for Eq. §5§ to hold we only re- 
quire that the distribution %j>\ can be written in the form 



■0i(i - t'\t') = t' 1 \(t/t')0{t - t'), which holds as long 
as there is no time scale in the problem other than the 
arrival time t' . Moreover, we require that G(p) can be 
Taylor expanded. An example when these two conditions 
are satisfied is the random walk model for transitions be- 
tween energy minima in a simple glass proposed by Ange- 



lani et al. [23j . These authors found that the rate of tran- 



sitions decays as c/t, where c > is a numerical constant. 
This decay rate of transitions corresponds to the waiting 
time distribution ipi(t~t'\t') = ct' c 9(t-t')/t 1+c , leading 
to the moment generating function G(p) = 1/(1 —p/c). 
Eqs. ([7]) to (fl3|) therefore apply to this process with 
fi = 1/c, g = 1/c 2 and n 3 — 2/c 3 , and thus also in 
this case the fluctuations vanish relatively to the mean. 

We finally consider the time average for a single real- 
ization n(t) of the counting process with a < 1. Such 
time averages are routinely measured in single particle 
tracking assays 14] • For the number of states this time 
average is defined in terms of the lag time A via 24] 



n(A) 



1 



h - A - 1 



ta-A 



1 Jti 



i(t + A) - n(t)]dt, (15) 



where the observation time interval of the trajectory is 
from t\ to t% ■ Averaging this behavior over many trajec- 
tories, the dominant behavior at ti,t% 3> to becomes 



n(A) 



A^hA 
ti — 1\ t\ 



(16) 



for A <t; ti — t\. This linear behavior in A contrasts the 
logarithmic time dependence of (n(t)). The linearity in A 
and the dependence on the measurement time is similar 
to the findings of weak ergodicity breaking of renewal 
CTRWs with scale- free waiting time density 14 , 2|| 2(| . 

We studied a general counting process in which suc- 
cessive transitions are affected by the age of the system 
through a local triggering process. This could be relevant 
for crack propagation in solids or counting statistics in en- 
sembles of systems with power-law transitions times such 
as blinking quantum dots. We presented analytical and 
numerical results for the temporal distribution to reach 
a state n as well as the fractional order moments of the 
process n(t). Interestingly the time evolution becomes 
logarithmic. Moreover, we found that the mean (n(t)} of 
the counting process grows faster than the deviations of 
n(t) for power-law waiting time densities with exponent 

< a < 1 and a > 2, while their ratio approaches a con- 
stant value in the intermediate case 1 < a < 2. We note 
that when the power-law exponent a > 2 the counting 
process becomes normal in the sense that the mean num- 
ber grows proportionally to the process time, (n(t)) ~ t, 
and in this respect corresponds to waiting time densities 
-0(t) with finite moments of any order, for instance, an 
exponential ip( T ) — ex P( — t/tq)/to. However, similar to 
biased subdiffusive CTRW [27|, the intermediate range 

1 < a < 2 still carries anomalous signatures. 
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Our results for the counting dynamics in an ageing 
system will be useful for the description of frustrated 
(e.g., glassy) systems, in which the renewal assumption 
of CTRW in an annealed environment no longer holds. 
As we showed, the generalized counting process exhibits 
weak ergodicity breaking. For the interpretation of single 
particle tracking experiments in such systems it will be 
interesting to obtain the full stochastic description of the 
position coordinate. We finally note that logarithmic sys- 
tem evolution was recently studied in the framework of 
superheavy-tailed waiting time densities within a renewal 
CTRW approach [28| . In contrast, our present results 
are based on a regular power-law waiting time density 
but with locally aged system update. 
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